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The positron-hydrogen system is considered, with the view to determining 
whether a scattering resonance below the positronium threshold is likely to 
occur, as has been suggested by Bransden and Jundi. The likeliest mechanism 
for producing such a resonance would be the existence of a virtual bound state 
in the positronium-proton channel, brought about by the large dipole polariz- 
ability and mass of the positronium atom. A modification of Hol^ien's technique 
enables one to reject spurious solutions without the use of the projection oper- 
ator technique. A non-variational adiabatic calculation seems to indicate that 
a resonance does exist, lying at 0.1 eV below the threshold. When, however, a 
variational treatment is employed, certain essential non-auiabatic terms appear. 
The resultant reduction in attraction serves to eliminate the apparent resonance. 
It is concluded that the previously reported resonance is probably an artifact. 
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I. INTRODUCTION 

Bransden and Jundi 1 have recently examined the positron-hydrogen system 
above the threshold for positronium formation (6.8 eV). With the aid of the 
M-Matrix method 2 they have extrapolated their results below the threshold, and 
have presented evidence that a resonance occurs in the e + -H elastic channel. 
Their estimate of the position of this resonance corresponds to a state of the 
positronium-proton system, bound by about 0.02 eV (1.6 x 10" 3 Rydberg). The 
present work is intended to examine the evidence for the existence of such a 
resonance. 

The simplest method of locating such a resonance would be the projection 
operator technique, 3 which has been very satisfactory in finding e"-H resonances 
below degenerate thresholds. 4 This technique effectively decouples the open and 
closed channel parts of the scattering wave function, and reduces the problem of 
resonances to that of bound states. The most likely physical mechanism which 
could produce binding is the electric polarization of the positronium atom in its 
ground state by the proton. Thus a wave function including this polarizability 
will be the basis of the present work. 

If we construct an appropriate scattering trial function 0 including polar- 
izability, the projection operator prescription for computing the resonant ener- 
gies is as follows: 

(1) Construct a projection operator Q which removes the hydrogen ground- 
state component from 0 , and 

(2) Carry out a Rayleigh-Ritz variational energy calculation using Q 0 as 


the trial function. 
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Using the simple form Q = 1 - |H>< H|, where |H> represents the hydrogen 
ground state r), x is the positron coordinate, and the trial function has the 
form 


-e- 

* 

II 

✓"S 

* 

U 

"9- 

H ( r ) + $ ( X - r )> 

(1) 

one can write the closed-channel part of 4> as 


Q'A = sK*. *) - 4> H (r) J 

d 3 r * 4> n (r*) <p( r, x), 

(la) 

Q4> = 3>(x, r) - 4> h (r) J 

d 3 r ' 0 H (r*) 4>(x, r') . 

(lb) 


If, for example, <l>were a close-coupling expansion in excited states of the target, 
then Qi <S>. We wish, however, to write $ in terms of the positronium atom 
in the field of a proton, so the second term in Eq. (lb) does not necessarily 
vanish. 

Hol^ien 5 has used a simplification of the projection operator method in 
searching for electron-atom resonances, and we will here apply a form of his 
approach to the positron problem. To do this, we first neglect the second term 
in Eq. (lb), and use <l> itself for the energy calculation. Some of the resulting 
energies will correspond to functions $ for which the integral in Eq. (lb) is 
very small; i.e., for which the projection onto the hydrogen ground state is very 
small. We follow Hol^ien in considering such states to be genuine resonances, 
although the energies thus obtained have no necessary upper bound character. 
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II. THE SCATTERING WAVE FUNCTION 

Let us consider the motion of a positronium atom (Ps) in the field of a 
fixed proton (P). We will use a trial function of the following form 

[x(R) + F (R) G (R , C . cos6>)] 0 (C), ( 2 ) 

where 0 (C) = (877) ~ 1/2 e _1/2 ^ is the positronium ground-state wave function, 
while R = 1/2 [x+r],{ = x- r and cos 9 - (R* £ )/RC- The functions x 
and F will be determined later, while G is to represent the adiabatic distor- 
tion of the positronium atom. 

The Hamiltonian can be written in atomic units, with energies in Rydbergs 
as 

H = H£ + V-1 -V 2 , (3) 

where 

H J =- 2 ( V ! + |)' [ H I + •§•]*«)=«• 


V = 2 






= v ( r * 0 (cos 9 ). 

0 ( odd ) 


The first-order (in V) adiabatic form for G is obtained by solving the equation 6 


[G, H t ] 0 = V0, 


( 4 ) 


and since by symmetry only odd multipoles are retained in the expansion for 
the potential, it is a good approximation to keep only the dipole [0=1] term. 
Then, G = g(R,C)P x (cos 8 ) and Eq. (4) becomes 


d 2 g 

dC 2 


+ 



dg 

dC 



( 5 ) 



where f = -£/R 2 for £ - 2R and f = -8R/ i 2 for £ - 2R. The properly con- 
tinuous, regular solution of Eq. (5) with continuous first derivative is found to 
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81 « 5 2R) = iF t £ + H -r ( 1 *s 1 [§ ' ’ (‘ f p )] 
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III. THE ADIABATIC NON- VARIATIONAL METHOD 

Our first evaluation of x and F will be non-variational, and resembles the 
polarized-orbital method 7 used previously for positron-hydrogen scattering. 8 
We assume that F = x in Eq. (2) and require x to satisfy the equation 

<([H - E] [1 + G]> X = 0- (7) 

where E = -1/2 + e , and the bracket is defined as 

<q(R.{)> =Jd 3 C 0(0 q(R. 5)^(0- (8) 

A bound state in the Ps - P channel would occur if a normalizable solution of 
Eq. (7) can be found with e < 0, and it would represent a resonance in the 
elastic e + -H channel near a positron energy 1/2 + e. Since 0, Eq. (7) 

can be written explicitly as 

[-y v l+ v 2 ] X = ^X. (9) 

where V 2 = <pV> is the adiabatic dipole potential, listed in Table I. Before 
solving the eigenvalue problem of Eq. (9), one can check whether any bound 
states exist by solving the S-Wave, zero-energy [e = 0] scattering problem. 
The number of nodes in the radial wave function will give the number of bound 
states. This radial function U = Ry is obtained by numerically integrating the 
equation 

- — U* + [V, - e] U = 0 (9a) 

2 2 

outward from the origin in the usual way. The result for e = 0 is plotted in 
Fig. 1, and one sees that two bound states occur. The more compact of these 
is spurious, and results from the omission of the second term in Eq. (lb), but 
the state giving the second node may represent the resonance found in Ref. (1). 
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The eigenvalues are'then obtained numerically from Eq. (9a), by integrating 
outward from the origin and inward from the asymptotic region. Continuity of 
the logarithmic derivative determines e . The lower eigenvalue is € ~-0.6 Ry; 
this energy is close to that of the IS state of hydrogen, and the corresponding 
eigenfunction peaks at R% 1.5, so this clearly is one of the spurious states 
which would have been eliminated by the Q operator technique. Thhvsecond 
eigenvalue is e = -7:5 x 10 " 3 Ry, and represents a resonance lying 0.1 eV be- 
low the positronium threshold. 9 In Fig. 2 the low-lying e + -e~ -P states are 
shown on a level diagram in order to clarify the situation. The two .bound 
states just discussed are in the closed channel below the ground state of posi- 
tronium [Ps(ls) + P]. The spurious state is seen to be in the vicinity of the 
ground state of hydrogen [H (Is) + e + ] , while the second state corresponds to 
a resonance in the open channel. In Section V we will discuss the relation be- 


tween this resonance and that found in Ref. 1. 
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IV. THE VARIATIONAL METHOD 

The assumptions made in Section III [that F = y and that Eq. (7) holds] de- 
pend on the ability of the positronium atom to adjust adiabatically to the polar- 
ization induced by the electric field of the proton. 10 This in turn depends 
roughly on the smallness of the velocity of the positronium center of mass. 
The very tightly bound (spurious) lower state probably has too much kinetic 
energy to satisfy the adiabatic conditions, but the resonance state may satisfy 
them. 


The variational method relaxes the above condition on F and determines 
both functions by free variation. Eq. (7) is replaced by the following two equa- 
tions: 


<[H -E] [x +FG]> = 0 


<G [H -E] tx +FG]> = 0, 


(10) 


whose solution corresponds to solving the Schrodinger equation in the restricted 
subspace spanned by the orthogonal vectors 4 > and G<£. Explicitly, these result 
in the following system of differential equations: 11 


[y V R + e ] x - V = 0 

[n U v H 


y 2 , e \ . y _ vv - V — 
V R + t J + V 2 " N dR 


f-v 2 x = o, 


where 


N= <G 2 > W--1<G V2 G >, V 2 -<GV> 


( 11 ) 


[These fimctions were evaluated numerically and are given in Table I.] Except 



for R - 00 , the relation F = y does not satisfy Eq. (11). Interesting features of 
Eq. (11) are the short-range potential W and the velocity-dependent term in- 
volving v N . 
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Letting y = R 1 U and F = R' 1 g one can write the S-wave , e = 0 form of 
Eq. (11) as follows: 


i-U" = V 2 g 


^-g' -V N g' + lv 2 -W t V„/R] g = V 2 u. 
£ 


(12) 


These equations resemble those of a two-state close-coupling calculation, 
and to treat them similar techniques are used. 12 Two linearly independent 
regular solutions of Eq. (12) exist for small R; their leading terms are 


and 


Uj =R, gl = 0 
U 2 = 0, g 2 = R • 


(13) 


Each of these solutions is integrated numerically from the origin out to a con- 
venient intermediate point R , at which point a linear combination 


U< (R) = AU X (R) + BU 2 (R) (14) 

with undetermined coefficients represents the numerical value of U. A similar 
combination is formed for g. Three linearly independent regular solutions of 
Eq. (12) exist for la^ge R. To find them, we use the asymptotic forms of the 
potentials given in Table I. Eq. (12) becomes 


1 _ 

2 


— U" = - 36 R~ 4 


~ g" - 2R _1 g‘-i! g = -J® U, 


1 

2 


43 


43 


(15) 


where terms of order R 2 were neglected in the second line. The three solu- 
tions required for large R are 
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U 3 = 1 - 6/R 2 , g 3 = U 3 , 

U 4 = R - 18/R , g 4 = U 4 - 43/9R , (16) 

U 5 = 0, g 5 = e-^ R £l + ySR + I(/3R) 2 j, (/ 3 2 = 36/43) . 

[Solutions 3 and 4 contain the first two terms of a series in inverse powers of R]. 
Each of these solutions is integrated numerically inward from an asymptotic 
point R Q (= 20 in the present work) to the intermediate point R where the linear 
combination 

U> (R) = U 3 (R) + CU 4 (R) + DU 5 (R) (17) 

represents the numerical value of U, and similarly for g. (The coefficient of 
U 3 is taken equal to unity, and U has the conventional normalization to unit 
slope at R - oo). Four equations represent the matching conditions at R = R: 

U< = U>, g< = g>, U'< - U> , g'< = g'>. (18) 

These determine the four constants A, B, C, D through a set of four linear 
algebraic equations. The function U(R) thus obtained is plotted on Fig. 1. This 
is the principal result of the present work: The node representing the spurious 
state occurs at larger R than for the non- variational case, and this decrease 
in effective attraction is sufficient to destroy the second node. 13 The variational 
treatment thus does not give any indication of a resonance below the Ps threshold. 
The radically different behavior produced by the inclusion of the non-adiabatic 
terms leads us to conclude that the resonance found in Ref. (1) is probably spurious. 
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V. DISCUSSION AND CONCLUSIONS 

To explain the relation between the present work and that of Ref. 1, let us 

recall that the latter considered the scattering of positrons on hydrogen above 

the threshold for positronium formation, using a polarized orbital type of wave 
+ 

function in both e -H and Ps-P channels. In our notation 


'P = Y] ( x ) tl + G„(x, r )^ ( r ) 


(i9) 


+ X 2 (R) [1 + G ps (R » 5 ) ] 0p s CO- 

If the Schrodinger equation [H-E]0 = 0 is projected onto c£ H (r) and <£ Ps (£) in 
turn, the coupled equations used in Ref. 1 are obtained [their Eq. (16)]. These 
are as follows: 


[V2 + k2 - V 1 (x) - V 2H (x)] Xj (x) = j d 3 r <£ H (r) 

[_V2 - k 3 + V (r, x)] [1 +G ps (R, C r ) *i ( x ) = 0 

[ V R+ k 2 - 2V 2P s ( R )^ 2 (R)= J d3 5 * P .<0 

[_V| -k* + 2V(R,£)] [1 +G h (x.r)] 0 Pl (Ox 2 (R) = O 


where the second-order potentials are the adiabatic polarization potentials of 
hydrogen [including several multipoles] and positronium [dipole distortion only], 
and where 


V (r , x) 


2 


_ 1 __ 1 _ 

x j x - r 


(21) 


and V (R.i )is given in Eq. (3). [A further neglect of all terms in G H , G Ps on the 
right side of Eq. (20) does not concern us here]. If Eq. (20) is now solved, the 
eigenphases or reaction matrix satisfy a minimum principle, but only if certain 
non-adiabatic terms are added. These are complicated and involve cross-channel 
couplings, but for the diagonal Ps term the form is 
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J”d 3 R X 2 («) [~| N (V| + k 2 2 ) + W + V N A ] X2 (R) , (22) 

where N, W, V N were defined in Eq. (11), and a similar correction term appears 
in the diagonal H term. The authors of Ref. 1 emphasize that, since they have 
omitted the non-adiabatic correction terms their results cannot be considered 
quantitative. We have here attempted, by including terms of this type in an 
optimal variational calculation, to indicate that the most likely result of the 
quantitative extension of the work of Ref. 1 would be the disappearance of the 
resonance. 

Recently, another above-threshold calculation has been reported, 14 which 
does not obtain the resonance discussed in Ref. (1). The method is very similar 
to that of Ref. (1), including both ground- state hydrogen and Ps as well as adi- 
abatic polarization potentials. A unique coordinate system, introduced into the 
Ps-P channel to simplify the form of the coupled equations, also provides a 
spurious centrifugal potential. It seems that this extra potential serves inad- 
vertently to simulate the effect of the non-adiabatic terms we have discussed 
above, and hence to eliminate the resonance. 

An interesting question may be raised now concerning the Ps - He + channel 
in e + -He scattering. This system is the same as Ps - P at large distances, but 
its additional Van der Waals attraction might be just enough to produce the 
resonance. The cross-section near threshold for e + -He scattering has been 
measured, 15 although the analysis of the experiment is not straightforward. It 
is possible that a resonance near threshold is needed to bring about agreement 
between experiment 15 and theory. 16 
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TABLE I 

Potentials Needed for the Numerical Solution of the Scattering Problem 


R 

-v 2 

N 

v N 

W 

0 

0 

0 

0 

0 

0.2 

0.282 

0.198 

-0.920 

0.801 

0.4 

0.745 

0.661 

-1.300 

1.944 

0.6 

1.108 

1.173 

-1.197 

2.613 

0.8 

1.306 

1.586 

-0.845 

2.752 

1.0 

1.358 

1.842 

-0.436 

2.536 

1.5 

1.135 

1.873 

+0.276 

1.521 

2.0 

0.785 

1.469 

0.464 

0.729 

2.5 

0.504 

1.026 

0.403 

0.314 

3.0 

0.316 

0.680 

0.288 

0.128 

4.0 

0.127 

0.291 

0.121 

1.97 (-2) 

5.0 

5.60 (-2) 

0.132 

4.92 (-2) 

2.90 (-3) 

7.5 

1.14 (-2) 

2.71 (-2) 

7.21 (-3) 

2.32 (-5) 

10.0 

3.60 (-3) 

8.60 (-3) 

1.72 (-3) 

1.76 (-7) 

R-*oo 

-36 R -4 

-86 R- 4 

-172 R~ s 

Exponential 
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FIGURE CAPTIONS 

Fig. 1. Radial wave- Functions U= Rx, for e = 0. Curve A is non-variational 
and shows two- nodes, while curve B is variational and has only one. 

The functions are normalized to unit slope asymptotically, with R 
measured in Bohr radii. 

Fig. 2. Energy level diagram for the e + -e -P system. The two bound states 

found in the non-variational method are shown: A is the spurious level, 
and B respresents the resonance found below the Ps threshold. The 
diagonal hatching indicates the scattering continua. The terminology 
"open" or "closed" channel is meaningful only below E = -1/2. 







